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title pl i$div pk Long | 
ident /1-002/ : Edit wHM1002 
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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS U LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
ooRPORAT i on NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


eaeneneeeee eee eneeeaeannnnen 


(ARBRARRARSAASAALASARAAALAREAA EASA ERAS ARERR A SARE SERS R SASSER ESSA ASAE READS SS 


ePe Se Ge Ge Se Se Se Ge Ge Se Ge Ge Se Se Se Se Se Ge Se Ge Se 


+ 
+ 


SOOO CCOOSOOOOOOOOOOOOOOOOOOoOO 
SOSCOSCSCOSOSOSOSOSOSOSOSOSOSOSOOOSOOOOoOoOSoSo 
SOoOOoOSCOSCOSOSOSOSOOSOSOSOOSOOSOOOSOOSOOSoOO 


oo 
Ss 
SOOCOCCOOCOOSOSCOSOSSOSOSOSOSOSSSSSSSSOSSOSOSCSOSOSOCOSOOOCOOOOOOOOOOOOOOOOO 


routine: 
PLIS$DIV_PK_LONG 


ooo 


ooocoo 
Soooooooeo 


facility: 
VAX/VMS PL1 runtime Library. 
abstract: 
The routine is called when precision and scale Fegutrenents for 
the quotient imply ay Site precision division. 
only called when such mu 
when the divisor has a precision of. digit F 
(Call pli$div_pkshort if multiple precision division is 
required and the divisor has precision less than 30 decimal digits). 
author: Peter Baum 30-jun-1980 | 


modifications: 


tipt ne Me dy oP Ba de 
e recision IV1IS1ON 18 required an 
Felon 36 or 31 decimal a? i 


Runtime routine performs fixed decimal (packed decimal) division. | 
| 


ts. 


1-002 Bill Matthews 29-September-1982 | 
Invoke macros $defdat and rtshare instead of $defopr and share. 
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documentation file: Cpl1.doc.codegen] THEORY .MEM 
functional description: 
This routine calculates: 


z=x/y 

let a = scale(z) + scale(y) = scale(x) = 31 + greece? 
b = scale(z) + scale(y) = scale(x) + prec(x 
c = 31 = prec(x) 
d = 31 = prec(y) 


this routine is called if b > 31 and d < 2 
Prior aa or call: 


f c not 0 then shift x left by c. 
Thus x is a 31 digit packed decimal. 


input: 


O(ap) # of arguments. < 

4(ap) address of dividend (shifted left by c) 

8(ap) address of divisor ‘ 

12(ap) precision of divisor (high order bytes zeroed) 

16(ap) address of quotient 

20(ap) precision of quotient (high order bytes zeroed) 

24(ap) a as defined above(high order bytes zeroed) 
output: 


OOO OOO 000009 09.09 09 09 09 09 09 SII NI NINN NIAAA AAAAAOOOU 
MEWS OONOAUES WN 0 ODNOAUE WN OOONOUSWh—O00 
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quotient returned at address specified by 16(ap) 
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96 
97 
38 optimization notes: 
100 1) Optimized for speed, not space. 
0 13) ¢ Optimized for y > 0. 
0 1 ¢ ) Assumes speed for register to register operations are the same 
0 10 for byte operations and longword operations, 
0 104 4) Many packed instruction sequences were timed. Do not shanes 
0 193 unless actual tests are made to determine relative speed. 
136 Tests were made on 11/780 and Comet. 
0 1 
00 109 | ie as 
0 19 possible optimizations: 
4 16 1) currently we always calculate the next 15 digits each 
11 iteration and then truncate the last iteration as part of 
0 114 the final step. We might be able to go through making 
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a a ed ed od td ot od = oo 8 2 3) 2 8 2 = 8 = = 2 a 2 a a 
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16-SEP-1984 02:22:23 
att 96:98:63 
calculations with fewer digits on this last pass. 


2) the classical trade-off involving flow paths which 
could be merged (sometimes apis in things Like 
adding 0) or merged with switches or left alone (resulting 
in space-speed tradeoff in favor of space). Since we are 
running on_a virtual machine, even the speed factor is not 
obvious. This came up when remainder was 0, when borrow 
was not needed, when less than 15 digits where required 
for the last part of the quotient, and when the new 
algorithm had found precisely what the next 15 orgie 
were but did not have a value yet for the remainder (the 
eer ~ papcnate always have a remainder when the quotient 
s known). 


3) New algorithm - special case y2=0 


4) New ss a - optimize branches according to probable 
sign of RCH). 


variable use: 


size 
in 
variable digits use 

yl 15 High order digits of divisor. 

y2 16 Low order digits of divisor. 

a 31 Initially dividend, thereafter 
remainders of successive divide 
operations. 

2 gy ap) Quotient. 

z2 Temporarily holds trial low 
order digits of quotient. : 

t1 31 High order digits of the remainder. 

t2 31 Holds the 15 low order digits of the 
46 digit remainder. 31 digits for 

ossible later changes. 

t3 16 olds the low order digits of the 

a 31 remainder. 


not overlap their operands. 


register usage: 


register use 
r6 a = additional gigits of precision required 
r7 stky(sp) which holds a copy of divisor 


AX/VMS Macro Vv04-00 
PLIRTL.SRCIPLIPKDIVL.MAR; 1 


compar ers used because packed instructions can 
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-SEP=-1984 3 VAX/VMS e cro V04-00 Page 4 | 
-SEP=1 198 4: $6; 83 PLIRTL.S RCIP LIPKDIVL.MAR;1 " (1) > 
172 ; r8 (ap) recision of y 
009 178 ; r9 rs mse of additional digits of the quotient 
88 1% : 10 that are to be found for next step 
: r 
0900 128 : rij patna) © = precision of quotient 
4 177 ; 
000 178; 
500 179 ;-- 
000 180 ; 
0000 13) 3; stack offsets for work area 
0000 1 § : 
464 13 Soffset 0,,<- H 
0000 184 <,16>,- ix, 31 digits 
0000 185 <Stky{,8>,- tyl 15 digits 
0000 186 <stky2,9>,- sy2 16 digits 
0000 187 <stkz2,16>,- :z2 31 digits 
0000 188 <stkt1,16>,- 3t1 31 digits 
0000 189 <stkt2,16>,- ;t2 31 digits (15 digits used) 
0000 190 <stky, 16>,- ty. 31 digits 
0000 191 <stkt3,9>,- :t3 16 digits 
0000 136 <stkt4,16>,- :t4 31 digits 
0000 19 <stksign 1>,- isign of quotient, 2 bits 
0000 194 <stklen,0>,- ; Length of work area 
0000 195 > : 
0000 196; 
0000 197 ; parameter offsets 
0000 198 ; 
0000 199 $Soffset 4,,<- 3 
0000 200 <x>,° sx = aieieent by reference 
0000 201 <y>,> zy = divisor 7? reference 
0v00 $08 <py>,- leant yl by value 
0000 20 <2>,° iz = quet ent by reference 
0000 204 <p2z>,- sprec(z) by value 
0000 205 <consta>,- 7a as defined above 
0000 206 > i 
0000 207 ; 
0000 208 ; constant data area 
st $18 ; (plidata.mar not used because they cause longword displacements to be used) 
it 211 rtshare 
1¢ 0000 $16 one: packed +1 : 
OC 0001 2135 zero: acked +0 
000 14 ;**warningtt the following two data definitions must be contiguous 
98 i 15 nines: yte 9 i *#16 - : (must be followed by 10**15-1) 
99 99 99 99 000 16 ninel5: :packed +999999999999999 *1 
90 38 4 88 Ooo 217 bignine: .packed $$000000000000009999999999999999_; 10##16 - 1 
00 00 00 01 Baie 18 teni5: .packed +1000000000000000 310*#15 
99 99 99 09 8 ; 4 neg9: .packed -9999999999999999 77(108##16 = 1) 
8 4 1 ; local symbol definitions 
O000000F 00 4 bytes_ to_sign=15 sbytes to sign for fixed decimal 31 
09 4 ; ; run time routine pli$div_pk_long 
CFFC $050 : z -entry pli$div_pk long, *M<iv,dv.r2.,r3.r4,r5,r6,r7,78,r9,r10,r11> 


PLISDIV_PK_LONG 16-SEP-1984 :22:23 VAX/VMS Macro v04-00 Page - 
S002 aes 96:96:67 PLIRTL.SRCIJPLIPKDIVL.MAR; 1 (1) 
5E 85 AE 9 F 8 movab -stklen(sp),s ;make room for temporaries 
28 18 AC Hi \ ; 9 mov consta cap) ig iget value of a 5 
7 51 AE 3 0 0 movab stky(sp),r? taddress of copy of divisor 
58 OC AC OD 038 1 mov py(a ) 66 zprecision of divisor 
5A 10 AC 00 OO035F ¢ mov | z(ap),r1 save address of quotient 
5B 14 AC 00 p3 movl pz(aps,ri1 precision of quotient 
7A AE 3 be 4 clrb atte len (se? :clear sign flag 
6— 04 BC IF 4 4A 5 movp #31,ax(ap), (sp) smove x, set cond. code 
24 «14 «O04F § bgtr 50$ sbranch if x > 0 
1¢ 8619 0051 bliss 0 sbranch if x <0 
005 38 ; 
edie 
08 BC 58 AAAF 00 37 005 41° cmpp4 #0,zero,r8,ay(ap) ;divisor zero? | 
09 13 005A 42 beat oy ibranch if divide by 0 
5B 00 AO AF O00 00 F8 OQOSC st ashp #0,#0,zero,#0,r11,(r10) ;z = 0 
04 0064 44 ret 
95 AF 00 98 AF 7 27 oe 245 308: divp #0,zero,#0,zero,r11,(r10) ;cause divide by zero 
04 O06E ses ret 
OO6F 47 ; 
OO6F 248 ;x not 0 
OO6F 249 ; 
006F 250 40s: : 
7A AE 96 QO6F 251 incb stksign(sp) set low order bit 
OF AE 97 Bore 526 decb bytes_to_sign(sp) zx < 0 so make it positive 
0075 254 ;determine sign of y 
0075 255 :y may be 0 at this point 
0075 256 ;optimized for y>0 
0075 257 ; 
0075 258 50s: ; | 
67 08 BC 58 34 0075 £259 movp r8,ay(ap),(r7) smove y into temporary 
OD 18 QO7A 260 bgeq 60$ ; sbranch if so 
7A AE 96 QO7C 261 incb Spkstontsp? ;set neg indicator | 
FF7A CF OV 08 BC 25 23 tha 262 subp6 = r 8, ay ap) ,#0,zero,r8, (r75 sconvert to positive 
0089 263; 
0089 264 ;get yl and y2 
0089 265 ; 
0089 266 60S: : | 
OF 00 67 58 + 4 F8 4:4 267 ashp #-16,78,(r7) ,#0,415,stky1 (sp) shigh order 15 digits of y 
1F 00 10 AE OF W F8 44 268 ashp #16, #15, stkyl(sp) ,#40,431,stktl(sp) sy with 16 low order zeroed 
. fF 6B Re fF & sie 269 subp6 #31,stkt1(sp),r8,(r7),#16,stky2(sp) ;16 low order digits y 
18 AE 00A 
00A4 70 ; 
BRAG 4 sprec(y) is large enough for yi to possibly not be 0 
10 AE OF FFS8 CF O00 37 ae 78 capes #0,zero,#15,stky1(sp) — ; : 
03 13 OOAC 74 beq 80$ sbranch if yl is zero 
0098 ==31 Boer o? brw 200$ zbranch if yl is not zero 
081 7% 3: y2(16) holds all of y 
0081 78 ; 
00B1 279 80s: ; 


—E 9 
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5 tO0e ety 7 09598 563 PLIRTL.S RCIP PLIPKDIVL.MAR;1 . (1) 
| 
18 AE 10 6€ 1F 37 Q0B1 0 cmpp4 #31,(sp),#16,stky2(sp) ix< A 
13 19 $oB i gnpes a+ (se Tee tbranch ff oy Shift x by 15 is of | 
6A SB 6F IF I8AE 1 g 9089 é divp #16,st ky3(sp), #31, + (Sp) ritcriQ) sz=x/ 
7, DR WwW ae 5 peel 8 mulp #16,s thy (sp),r11,(r10 sa3tstktlisp) ptl=(x/y2) *y2 
ig 11. OOCA gee brb 1108 H 
5B OO FF2F CF 00 oy F8 Bate 85 95$: ashp #0,#0,zero,#0,r11,(r10) sclear quotient 
\é 11 0005 286 brb 115$ 
7~T He & ne 25 std 587 100$: mulp #16,stky2(sp) ,415,stkz2(sp) ,A#31,stktl(sp) ;tl=y2*z2 
6E 1F 31 AE 1F 0E1 88 110$: subp4 #31,stkt1(sp) 431, (sp) exe =x-t1 
3E 3 GOE? 59 beql 130 . Pr ibranch if remainder = 0 
it 4 sdetermine r, the number of the next low order digits to obtain | 
59 OF 00 O0E9 38 115$: movi #15,r9 pr=15 
59 56 01 OOEC 94 cmpl r6,r :a>15? 
03 18 OOEF 95 bgeq 1308, anes if larger 
59 56 DO OOF1 296 movl reer 
31 AE 1F OO 6€ 1F 59 £F8 OOF& 297 130$: ashp r31, (sp), be od sstktt sp) jo x left by r 
31 AE 1F «6934 «(OOFC «6298 movp a3, stkel (sp) py ba ack into x 
31 AE 5B 00 6A 5B 59 FB 0101 299 ashp rit, (r10 #0 ae sstkel(Sp} sshift z left by r 
6A 31 AE 5B 34 0109 300 movp aty stktl(sp),(r ack in 
OF 6& 1F 18 ne 19 27 0108 01 divp #16, *Stky2(sp), #31, (sp), wis, otk 22 (sp) 122015)=x/y2 
6A 5B 21 AE OF 20 Q117 308 addp4 s- #15, stkz2(sp),r11,(r10) ;2=z+22 | 
56 59 C2 Q11D 30 subl2 r9,r6 azar 
BS 12 0120 304 bneq 100$ sbranch if more 
13 7A AE EB Q122 305 140$: blbs stksign(sp),155$ tbranch if quotient <0 
04 0126 306 ret : 
0127 307 ; , 
0127 308 ;remainder = 0 | 
0127 309 ; 
0127 310 150s: sremainder = 0 
31 AE 465B) O00 SGA CSB.) SOS6SCOFBCO27—s 311 ashp —s 6,711, (r10),#0,r11 sstktt (sp) saccount for scale 
OB 7A AE €E8 bigs 31§ blbs stksign(sp) ‘3 ranch if quotient is negative 
6A 31 AE 58 be O13 $12 aove rii,stktl(sp),(r10) ; Copy back into quotient 
re 
31 AE 6A 5B 34 Q139 +315 155$: movp r11,(r10),stktl(sp) copy quotient into temporary 
O13E 316 160$: nter if ti holds quotient 
5B FEBB CF O00 31 AE 7 23 b125 317 subp6 rill ,stktl(sp),40,zero, ric P10) smake z negative 
04 014 1 ret 
gia8 15 ; | 
0149 20 ;New division algorithm 
0149 $! . 
0149 é 2008: : 
0149 § : 
Spe : sinsure that we have x<y 
58 6€ 1F 37 Q149 6 cmpp4 #31,(sp),r8,(r7) apey? 
19 19 biaF 7 biss e reneh iM x<y 
1F 667)~—(58 32 15 8 divp (r7) #31 (sp) rit (rb). :2=x/ 
6A 5B 25 0157 9 mulp rif, (16) 75 #31, stktlisp); tM: =ytz | 
31 AE TF <2 15F 0 subp4 #31 stkttlops. 51 erie :x=xet 
cO 3613) «(0165 1 beq 1508 sbranch if remainder is zero 
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#0,#0,zero,#0,r11,(r10) :clear quotient 


i Assumes z is defined 
<x< 


and ré=a gives the number of additional digits required. 
:Determine r = # of digits for next part of quotient. 


vw OW 
oo "0 


tcalculate 22 


calculate yeree 
t3 = 
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r6,r9 irza 

r6,#15 ta>15? 

246$ tbranch if larger 
#15,°9 tre] 

r9,r6 jupdate ré6 = 


= min(B-1,Cx/y1)) 


#15, stky1(sp) #31, (sp), #31 »stktl(sp);tl=x/y1 
$7) stktl(sp) ,#16,nines et1>108#16-12 
ranch if not 


310 
#15,ninel5,stkz2(sp) ;move in base-1 
320$ iskip ashp, you have 15 digits 


high order 16 digits of y2*z2 

low order 15 digits of y2*z2 

#-1,431,stkt1l(sp) ,40,415,stkz2(sp) swe only get 15 digits 
#15,stkz2(sp) ,#16,stky2(sp) #31,stktl(sp) ;tl=y2*z2 

#-15, #31, stktl(sp) #0,#16,stkt3(sp) ;t3(16)=t1 shifted right 15 
#15, #16,stkt3(sp) #0,431,stkt4(sp) ;t4(31) = t3 shifted left 15 
#31,stkt4(sp) ,#31,stktl(sp) A#15,stkt2(sp) ;t2(15)=t1-t4 

330$ sbranch if no borrow required 


a, 


“31 high order digits of x(46) - y*z2 
5 low order digits of x(46) = y*z2 


#15, stkz2(sp) #15, stkyl(sp) #30, stkt4(sp) ;t4(30)=y1*z2 
#1,#30,stkt4(sp) ,#0,431,stkt1(sp);31 high order of 46 


#16,stkt3(sp) 431, a 3t1(31)= eatetoesgcte) 
#31,stktl(sp), #31. (sp) ,# stkt4(sp) 7t4=x-tl 


325$ sbranch if R(H) negative (t4 < = 0) 
¥i,one,#31,stkt4(sp) ,.#31.stkti(sp) ;tl=t4-1 (borrow 1) 
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OF FEIG CF 10 OF 23 0208 378 subp6 #15, stkt2(sp),#16,ten15,#15,stkt4(sp) ;t4=10*#15-t2 | 
41 AE 34 A 379 mov #15,stkt4(sp),stkt2(sp) scopy back into t2 
1 0910 $80 bro, 3708 ‘ vs ape 
31 AE 1F ef 1 4 325$:  movp #31,stkt4(sp),stktl(sp) ;make t1_ hold high order R(H) 
= Ht § brw 500$ zR(H) < 0 (can not be zero) | 
6 1B 54 ino borrow, t2 = 0 i 
1B 82 zcalculate R(H) = t1(31) = 31 “ie order digits of x(46) = y*z2 
0 18 86 ; t2(15) = 0 = 15 low order digits of x(46) = y*z2 | 
0218 88 $308: 
25 8 18 9 mulp #15,stkz2(sp) ,415,stkyl(sp) ,#30,stktl(sp) ;t1(30)=22(15) *y1(15) | 
F8 0295 390 ashp #1,430,stktl(sp) ,40,431,stkt4(sp) ;31 high order 46 
2 O93 91 addp4 #16.stkt3(sp) 431 etkes (sp) 3t4=t3(16)4+t4 | 
2 OS 92 subp6 #31,stkt4(sp),431,(sp) ,#3i,stktl(sp) ;tl=x-t4 
12 033 sar bneq 370$ sbranch if remainder not zero | 
0240 95 Iti titi iii i ttt i iii iii iti iii iii 
0240 396 ;:« * | 
0240 397 ;* R(H) = 0, i.e. * 
0240 398 :*remainder is zero ’ * 
0240 399 ;*2z2(15) holds last non-zero digits of the quotient * 
bse8 ret :*quotient has not been shifted yet to recieve 22 * 
4 ® 
0240 402 ;* this is a local routine with no exit other than ret * 
0240 403 peeee eee eee eee eee EERE AAA AREER AAA A AAA AAA AAA AEE Eee 
0240 404; 
0240 405 $408: 
F8 0240 406 ashp r9,r11,(r10) ,#0.r11,stkt4(sp) ;t4=z shifted left by r9=r 
34 0248 407 mov rii,stkt4(sp), (r10) COPY back quotient E 
C2 0240 408 subl2 .f sshift needed to leave r9 digits 
F8 6520 409 ashp r9,#15,stkz2(sp) ,#0,415,stktl(sp); low order digits of z | 
20 0259 410 addp4 #15, stkt1(sp),r11,(r10) ;z=2+22 
DS O25F 411 tst r6 ja= 0? , 
13 0261 tig beql 395$ branch if a=0 (last iteration) 
F8 0263 41 ashp r6,r11,(r10) ,40 ril,stkti(sp);ad ust for scale 
E8 0268 414 350$: blbs  stksign(sp),410$ :branch if quotient < 0 
34 OQO26F 415 movp rii,stktl(sp),(r10) sback into quotient | 
04 Bs ee $18 ret : 
19 3 rt 370$: blss 500$ sbranch if R(H)<0 
8 77 419 Seeeee eee eee eee CeCe T AeA AT AAA A AeA AAA AAA AAA TATE | 
77 «420 :*R(H) > 0, ie. * | 
0277 «421 3;* t1>0 * | 
0277 4 ¢ z% * 
0 77 4 PRRARA AAR AREA AAA AAA AAA A AAA AAA AAA ARATE AAA ARATE A AAR AA AAA A RAAT TRATES 
77) 0 426 
F8 77 425 $80$: ashp r9,r11,(r10) #O.rll,stkt4(sp) ;t4=z shifted left by r9=r 
34 7F OG § mov rii,stkt4(sps, (r10) copy back quotient | 
0 B46 4 tst c 3a =0?. : 
1 86 «6428 beql 390$ sbranch if a=0 (last iteration) 
0288 429; 


PLISDIV_PK_LONG 
1-002 
6A 8 21 AE OF 
1f 6000 (31 AE. SOF OOF 
6A AE 
1F 441 AE OOF «6A AE IF 
6£ 
FECF 
59 F 
OF 00 21 AE OF 59 
31 AE 
6A SB 31 AE OF 
01 7A AE 
31 AE 6A 5B 
5B 46FDSA CF «600006031 AE 5B 
6A 
FDS2 CF =610—Ss 31 AE IF 
1F 600) «631 AE OTF OOF 
6A AE 
1F 6A AE 1F 41 AE OF 
31 A 
na ow oe «(6S 
48 
3E 
21 AE OF FDOA CF 01 
6A AE 5B 00 6A SB 59 
6A AE 58 
56 
OE 
21 AE OF 
1 AE 1F 
FESA 
59 of 
AE OF 59 


ny Tw 
Ww coon 


OWN 
WMS @rwwoo 


WWD 
—2£O 


PAA OPE EE EWA 


Foe oP ot ot ot ot ot ot ot ot et et a oe a a 
WONAULWN—OOON AUEWR— COWODNAUN 


460 


Foe oS oP ot oY oP ot ot ot ot 
ANNO 


SOOOOOCOCOCOOCOOSOOCOCOCOOOCSCOCOCOCOSOO OOOO OCOCOOCO COO OCOOCOOCOOCOOOOOOOOOOOOOO 
WONAULSWNM—O”O 


AA AAANAA NAAN ANI IPOPPIPYPYINYININININYNININPIDININIPIPPYIPYNYIPININIPIPIPIPOPINPINYPOPOPONPONPOPOPINYPONPTG 


i SS ss BQO DDOOOO BRBMMMMNIGIOVCIVANAAHOAHOOOOOCCWWWWIW Yr KY Ky KY KP Yr PP OOOMWMwo 
MRWOcODOWNOY bP POO OS TMONMVIO DLP KS Yr KwrKwrKr Kr Kr Kr rOO Ph Pr OUIMIOOWWIWWWNO ™ ~UIT OOOO 


RCH) <0 
3* so check if R( H) + Y > = 0 * 
® 


. 
FRR ER RAE EATERRE REE 


§00$: 


430 ;a not 0 
d 5} ddp4 
a 
4 : ash 
434 addpé 
: brw 
ta = 0 
$90$:  subl2 
ashp 
addp4 
395$: blbs 
ret 
400$: 
movp 
410$: subp6 
ret 


AX/VMS Macro VO 


Mor Seen18ke $9:99507 EPRINTE eae TOL YORDGOL.MAR; 1 


#15, ghze(so) o711,(r10) ;z=2+22 
W15,#351,stktl(sp), #0, #31 stkte (sp) it4=t1 shifted left 15 


#31,stkt4(sp) ,#15,stkt2(sp) ,#31,(sp) ;x=t4+t2 
230$ 


#15,r ishift needed to leave r9 sigits 
Hs #8, stkz2(sp),#0,415,stktl(sp);low order digits o 


#15,stktl(sp),r11,(r10) ;2z=2+22 
stksign(sp), 460s sbranch if quotient <0 


r11,(r10),stkt1(sp) 


onghe copy of quotient 
rit ‘stktl(sp), #0,zero.rli, 


(r10) ;make z negative 


2 Wiest it" tri peti chi 3 pn a ci ght ah aan ny eh Sa poke 


& 
* 


® 


#31,stktl(sp),#16,neg9 ;ti< Bee dy righ -1) . 


tbranch +Y <0 
#15,431,stkt1(sp) ,40, #31, stkt4 (sp) vanite 1 left 15 
#15 ,stkt2(sp) ,M31,stkt4(sp) #31,stktl(sp) ;tl=t1-t2 


r8,(r7) ,#31,stktl(sp) ptlstity 
600$ sbrench” if R(H) + ¥ <0 


530$ 
#1 sone. #15,stk 52(sp) 
r9.r11.(r10),#0,r1 


| an, vie z shifted left by r9=r 
rif aabasteas, (r10) 


scopy, back quotient 
7 RAL, if a=0 (last iteration) 


p).ril 


(r10) ;z=2+22 
5)! * (sp) 1 


;x=t 


9 shift needed to leave r9 digits 
5,stkz2(sp) ,#0,415 Athtilees: low order digits o 


:branch, if no remainder (first adjust 22) 


PLISDIV_PK_LONG 16-SEP-1984 0 723° ~VAX/ vas 8 cro V04-00 Page 
be 6-SEP pry ope 4 $8 °67 PLIRTL.S RCIP LIPKDIVL.MAR;1 ° 
6A 5B 31-AE_ OF 20 480 addp4 = #15, stkt1(sp),r11,(r10) ;z=z+22 
8C 7A AE 3 481 blbs tksign(sp) ,460$ sbranch if quotient <0 
4 ? ¥ ret 
4 4 i jrenetnder is zero 
21 AE «OOF:)=«C FCCC CRO 32 486 $308 subp4 #1,one,#15,stkz2(sp) 322=22-1 
FFO6 1 1 brw 346$ 3go add to quotient, then ret 
rts 5 Suneeenwenenseorsnenescenesenenecenteceseoqnncessneseceeseeeeoenenseess 
3* t 
491 s*R(H) + ¥ < 0 tr 
298 | s*calculate L = ee ie CxX/(y1+1)]) > 
187 - by theorem 5, L=(X/( yield * 
*® 
495 ;: PUTT T TTT TITEL TTT TTT TTT TTT TTL TT TTT TT TTT TTT 


497 600$: addp6 #1,one,#15,stkyl(sp),#16,stkt3(sp) ;t3=y1+1 
498 divp #16,stkt3(sp) 431, (sp) ,#16,stktl(sp) ;tl=x/t3 


10 10 AE OF  FCC1 CF : 01 21 
10 6€ <1 61AE_ 10 27 
31 A 


13 t3 = high order 16 digits of ygel 
3 t2 = low order 15 di 


99 ; 

00 isplit up y2et 

055.80 FS. 

03 ashp #-1,#16,stkti(sp), o. 15,stkz2(sp) ;we only get 15 digits 


OF OO 31 AE 10 8F FB 


Tt Be WO Ae Ff DB mulp #15,stkz2(sp) ,A#16,stky2(sp) .#31,stktl(sp) ;tl=y2*z2 

10 00 31 AE TF 4 F8 505 ashp #-15,H31,stktl(sp) #0, 416,stkt3(sp) ;t3(16)=t1 shifted right 15 
1F 600 661 AE 10 ‘ 4 F8 506 ashp #15, #16, stkt3(sp) #0,431,stkt4(sp) ;t4(31) = t3 shifted left 15 

OF 31 AE 1F 6A ne ps 23 507 subp6 «= #31,stkt4(sp),#31,stkt1l(sp) ,#15,stkt2(sp) ;t2(15)=t1-t4 


WO BUC B 2 SS SS SS SS HHH AH AMNAHNGUWNO PWOCHWAMMMMMNOUWP PPP PPP PP PNM HMyS & 
uw 
o 
7 


7 

7 

7 

7 508 ; 

f 203 FRR REAR ETE EERE EAE ARE EEE 
7* as 

4 ai): 3* calculate R(L) ~ 

4 21¢ eneneeeeennceesneeeneneeeessennenenstenernenneeqoneenseqnreenseennneees 

14; 
41 13 f 315 beql 630$ zbranch if no borrow required 

8 213 tborrow is -1, t2 not 0 

8 518 calculate R(L) = 

8 519 ; tat) = 31 high order digits of x(46) - y*z2 

8 520 ; t2015) = 15 low order qigits of of x(46) = y*z2 

8 521 ;note: it is always true that R(L 

OF FC92 CF 10 41 a e 23 3 5 . subp6 #15,stkt2(sp),#16,ten15,415, 5 het 7t4=10**15-t2 
41 AE ga AE OF 34 i 286 movp #15, stkt4(sp), stkt2(sp) scopy we 
1E 10 AE OF 1AE OF 25 3 525 mulp #15,stkz2(sp), #15 »stkyl(sp) ,#3 Ay 3t4(30)=y1*22 
9 
A 


AA 
1F O00 6A AE re. 01 = =F8 


“ 
~ 
o 


ashp #1, #30, stkt4(sp) ,#0,431,stktl(sp) ;high order 31 of 46 


Onmrmw non 
—~— 
> 
m 
CSOOOCOOCOCOCOCOCOCOCCOCOCOCOCOOCOS SOOO COOSCOCSCOCOOCCOCOOCOOOCOCOOCOOCOOOOOOCOOOO 


59 
PLISDIV_PK_LONG 16-SEP-1984 3 VAX/VMS Macro V04-00 p 11 
Sb02 -SEP=1984 tf $6: 63 LPLIRTLeSRESPLIPKDIVL.MAR;1 oo (1 


SNe ee 


31 AE IF _-—\. $9 AS 287 addp4 OUS-stht step} «2 jetkel (9p) ESF a toatl eri per sere) 
~~ & ww 1 < ph os 8 subp6 1istktl(sp) #31, (sp) ,#31,stkt4(sp) ;t4=x-tl 
1F 6A AE 1F FC4& OF 01 23 ; 529 subp6 «#1, one, #31,stkt4(sp) ,#31,stktl(sp) ;tl=t4-1 (borrow 1) 
: 
s 11 ‘ - brb 700$ 3 
C § ‘no borrow, t2 = 0 
: 7 ;calculate Ric) = t1(31) = x = yl*z2 
C 5 4308: 
1€ 10AE OF 21 ne 4 25 es 6 mulp #15,stkz2(sp) #15 ,stkyl(sp) #30, stktl(sp) ;t1(30)=y1*z2 
1F OO 31 AE ". a F8 8 " 537 ashp #1,#30,stkt1(sp),40,431,stkt4(sp) shigh order 31 of 46 
oo AE 1F 61AE 10 39 0305 538 addp4 - #16, stkt3(sp) ,4#31, ogthee tae) 304(31)=t4(31)+t3(16) 
6— 1F 6A AE_ IF 0O30C 3=— 539 subp6 #31, stkt4(sp), "#31. (sp) ,A#31,stktl(sp) ;tl=x-t4 
31 AE 03E3 
0O3E5 540; 
O3E5 541 PRRA AAA AAA AAA AAT A AA AAA ATTA AA AATATAAAAATAAA AAA AAA A AAA e eee eee 
O3E5 276 3* * 
O3E5 543 ;* calculate 22 = L + R(L)/Y * 
O3E5 544 ;* 4 
Oeee ot2 FRR EAR AAA EERE AAA RARE AERA EAE EEA E ARERR EEE EE 
0 
1F 6000s 31 AE A oF F8 O3E> 547 700$: ashp  #15,#31,stkt1(sp),#0,#31,stkt4(sp) ;shift t1 left 15 
6A AE 1F 41 AE OF 20 OSEE 548 addp4 = #15, stkt2(sp) #31, stkt4(sp) ;t4=ta+t 2 
OF 6A AE IF Oe zs 27 pare 549 divp r8,(r7),#31 stkt4(sp), #15, stktl(sp) ; t1=t4/y 
21 AE OF 31 AE OF $9 O3FE 550 addp4 #15, stkt1l(sp) .#15,stkz2(sp) ;z2=z2+t1 
FD9E 31 0405 551 brw 320$ ; 
0408 $25 
0408 £55 .end 


984 02:22:23 VAX/VMS Macro V04-00 Page 12 | 
3B 06:96:63 PLIRTL.SRCJPLIPKDIVL.MAR; 1 ° (1) 


K 9 
PLISDIV_PK_LONG 16-SEP=- 
6-SEP- 


1 
Symbol fFabTe 1 


| 
| 
BIGNINE 0000008 R 02 
BYTES TO.SIGN = Q0QQOOOF | 
CONSTA 0 BAe 
DIR.. = 00 1 
EGS 9000024 R 0 | 
NINE15 0 R 8 
NINES 00 R 
ONE 00 f 0 | 
PLISDIV_PK_LONG 0000002D RG 0 
PY 0000000C | 
00000014 
STKLEN B80 0076 | 
STKSIGN 000007A 
KT1 00000031 
STKT 00000041 
STKT 00000061 
STKTG 0000006A 
STKY 00000051 
STKY1 00000010 
STKY2 00000018 
STKZ2 00000021 
TEN15 00000018 R 02 
00000004 
Y 00000008 
Zz 00000010 
ZERO 00000001 R 02 
power eeen em eanne + 
H Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
. ° 00000000 ( 0.) 00 ¢ O.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS 00000078 ( 123.) 01 ¢ 1.) \NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
_PLISCODE 00000408 ( 1032.) 02 (¢ 2.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
poem oe mow ee me ee ee oo 
: Performance indicators : 
Phase Page faults CPU Time Elapsed Time 
Initialization 13 0:00:00.07 0:00:00.31 | 
Command processing 74 09:09:9 a 88 0:09:13 
Pass 1 92 00:00:02.36 Be BR: Bao 
Symbol table sort 0 00:00:00.05 0:00:00.05 
Pass 2 105 90:00:01.19 90 -00:02.95 
Symbol table output 2 88:8 :30.83 82038 39 
Psect synopsis output 1 00: 8: 0.0 0:00:00. 3 
Cross-reference output 9 38:8 320-8 808308 -8 
Assembler run totals 28 0:00:04. 0:00:10.09 


The working set Limit was 900 popes. 

15119 bytes (30 pages) of virtual memory were used to buffer the intermediate goge. 

There were 10 pages of symbol table space allocated to hold 26 non-local and 37 local symbols. 
553 source Lines were read in Pass 1, producing 14 object records in Pass 2. 


PLISDIV_PK_LONG 
VAX=11 Macro Run Statistics 


3 pages of virtual memory were used to define 3 macros. 


! Macro Library statistics : 


perenne ne ceo oe wee ew nem ewan $+ 


Macro library name Macros defined 


~$255$0UA28: PLIRTL sPerL IA Tnac MLB; 1 
55$DUA28 $¥SL183 STARLET. MLB; 2 

Torace (all Libraries) 

44 GETS were required to define 3 macros. 


There were no errors, warnings or information messages. 


AX/VMS Macro V04-00 P 
PLIRTL.SRCIJPLIPKDIVL.MAR; 1 


Pa19Be 11:89:01 


MACRO/ENABLE=SUPPRESSION/DISABLE=TRACEBACK/LIS=LIS$:PLIPKDIVL/OBJ=OBJ$:PLIPKDIVL MSRC$:PLIPKDIVL/UPDATE=(ENH$:PLIPKDIVL)+LIB$:PLIRTM 
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